Draft version January 15, 2013 

Preprint typeset using L£-T^X style cmulatcapj v. 11/26/03 



EXTREMELY BROAD RADIO RECOMBINATION MASER LINES TOWARD THE HIGH- VELOCITY 

IONIZED JET IN CEPHEUS A HW2 

I. Jimenez-Serra 1 , J. Marti'n-Pintado 2 , A. Baez-Rubio 2 , N. Patel 1 and C. Thum 3 

Draft version January 15, 2013 



o 

(N 



(N 

a,: 
6 

CO 



> 
o 

o 



ABSTRACT 

We present the first detection of the H40a, H34a and H31a radio recombination lines (RRLs) at 
millimeter wavelengths toward the high-velocity, ionized jet in the Cepheus A HW2 star forming 
region. From our single-dish and interfcromctric observations, we find that the measured RRLs show 
extremely broad asymmetric line profiles with zero-intensity linewidths of ~1100kms -1 . From the 
linewidths, we estimate a terminal velocity for the ionized gas in the jet of >500kms _1 , consistent 
with that obtained from the proper motions of the HW2 radio jet. The total integrated line-to- 
continuum flux ratios of the H40a, H34a and H31a lines are 43, 229 and 280 km s -1 , clearly deviating 
from LTE predictions. These ratios are very similar to those observed for the RRL maser toward 
MWC349A, suggesting that the intensities of the RRLs toward HW2 are affected by maser emission. 
Our radiative transfer modeling of the RRLs shows that their asymmetric profiles could be explained 
by maser emission arising from a bi-conical radio jet with a semi-aperture angle of 18°, electron density 
distribution varying as r and turbulent and expanding wind velocities of 60 and 500 km s . 
Subject headings: stars: formation — masers — ISM: individual (Cepheus A) — ISM: jets and outflows 



1. INTRODUCTION 

Radio recombination lines (RRLs) are excellent probes 
of the kinematics of the ionized gas in ultracompact 
(UC) HII regions (Garay 1989; Churchwell et al. 1989). 
These lines typically show simple Gaussian profiles with 
linewidths of 25-30 km s _1 , attributed to unresolved gas 
motions and/or pressure broadening (Gaume et al. 1995; 
Afflerbach et al. 1996; Keto et al. 2008). 

In addition to the classical UC HII regions with 
simple narrow (^30kms _1 ) Gaussian RRLs, Jaffe & 
Martin-Pintado (1999) reported that 30% of the ob- 
served UC HII regions show even broader RRL emission 
with linewidths of 70-200 km s -1 . These sources have 
power law continuum spectra with spectral indeces ~0.6 
(characteristic of constant velocity stellar winds; Olnon 
1975), and elongated/bipolar morphologies resembling 
ionized flows. Jaffe & Martin-Pintado (1999) proposed 
that broad RRL emission could arise from bipolar ion- 
ized winds generated in photo-evaporating neutral disks 
(Hollenbach et al. 1994; Gorti & Hollenbach 2009). 

MWC349A is the best studied object with broad H 
recombination lines at optical, IR, centimeter and mil- 
limeter wavelengths (Hartmann, Jaffe & Huchra 1980; 
Hamann & Simon 1986; Altenhoff, Strittmatter & Wend- 
ker 1981; M art in- Pint ado et al. 1989). This source, which 
has a rotating edge-on disk with a bipolar ionized flow 
(Cohen et al. 1985), is unique in its category because its 
RRLs at A<2 mm (quantum numbers n<35) are masers 
(optical depths <-l; Martin-Pintado et al. 1989, 1994). 
The RRL maser spots are located on the ionized sur- 
face of the disk (Planesas, Martin-Pintado & Serabyn 
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1992; Weintroub et al. 2008; Mart in- Pint ado et al. 2011), 
where the densities of the ionized gas are high enough 
(>10 6 cm~ 3 ) to invert the population of the levels in- 
volved in the RRL (Walmsley 1990). Although these 
masers are expected to be found in other UC HII re- 
gions (Martin-Pintado 2002), MWC349 is the only ob- 
ject where this emission has been reported so far in star 
forming regions. 

We present the first detection of RRL maser emis- 
sion toward the Cepheus A HW2 high-mass star form- 
ing region. This source shows a collimated, high- 
velocity ionized jet, with a continuum spectral index 
~0.6 (Rodriguez et al. 1994). Like MWC349A, the 
molecular material around HW2 is mainly distributed in 
a neutral circumstellar disk (Patel et al. 2005; Jimenez- 
Serra et al. 2009) that seems to be photo-evaporating 
(Jimenez-Serra et al. 2007). The proper motions of the 
HW2 radio jet suggest an expanding velocity for the out- 
flowing gas of >500kms _1 (Curiel et al. 2006). This is 
consistent with the extremely broad line profiles (zero- 
intensity linewidths of ^llOOkms -1 ) of the H40a, H34a 
and H31a maser lines detected in our single-dish and in- 
terfcrometric observations. Cepheus A HW2 is the sec- 
ond RRL maser object detected to date in star forming 
regions. 

2. OBSERVATIONS 

The H40a line at ^99 GHz was observed toward 
Cepheus A HW2 with the IRAM Plateau de Bure In- 
terferometer (PdBI) in the D configuration (beam of 
5.0"x4.5", P.A.=91 ) 4 . The phase center was set at 
a(J2000)=22 /t 56 m 17.98 s , S( J2000)=+62°01'49.5". The 
wide-band correlator WideX provided a continuous fre- 
quency coverage of 3.6 GHz, and a spectral resolution 
of 1.95 MHz (i.e. a velocity coverage and resolution of 

4 Based on observations carried out with the IRAM 30 m 
and Plateau de Bure Interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain) 
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Fig. 1. — Upper panels: Full spectra of the H40ct, H34a, and H31a RRLs measured toward Cepheus A HW2 with the PdBI, 30 m and 
SMA at 3, 1.9 and 1.4 mm, respectively. Box shows the region of the RRL spectra shown in the lower panels of this Figure. Dotted lines 
indicate the extension in velocity from vlsr^— 500 to 600 kms -1 of the H40o, H34«, and H31o RRLs. Lower panels: Zoom-in of the 
continuum-subtracted RRLs spectra, after being smoothed to a velocity resolution of ~15-17 kms -1 . Identifications of some molecular 
lines in the RRL spectra are shown. We also report the individual Gaussian fits to the blue-shifted and red-shifted velocity components of 
the RRL emission (dashed lines), and the total Gaussian fit to these lines (solid lines; see also Table 1). 



-10900 km s™ 1 and -6 kms" 1 , respectively). 3C454.3 
(23 Jy) and 3C273 (12 Jy) were used as bandpass cali- 
brators. We observed MWC349 (1.2 Jy) and 1749+096 
(2Jy) as flux calibrators, and 0016+731 (0.6 Jy) and 
2146+608 (0.3 Jy), as phase calibrators. Data reduction, 
calibration, imaging and cleaning were performed with 
the GILDAS software 5 . 

The H34a line at -160 GHz, was detected toward HW2 
with the IRAM 30 m telescope in a single-pointing ob- 
servation during 140 min. The EMIR El receiver was 
tuned to single sideband (SSB) with rejection of >10 dB. 
The beam size was 15". The wide-band auto-correlator 
WILMA provided a total bandwidth and spectral resolu- 
tion of 4 GHz and 2 MHz (-6700 kms" 1 and 3.7 kms™ 1 , 
respectively). Typical system temperatures were 200- 
225 K. Intensities were calibrated in Tj[, and converted 

5 See http://www.iram.fr/IRAMFR/GILDAS. 



into total flux (Jy) by using S/T^=6.4 Jy/K 6 . 

The H31a line at 210 GHz was imaged toward HW2 
with the Submillimcter Array (SMA) in the subcompact 
configuration (beam of 4.8"x3.8", P.A.=21 ) 7 . The cor- 
relator setup provided a total bandwidth of 4 GHz per 
side-band, and a spectral resolution of 0.8 MHz. This 
corresponds to a velocity coverage of —5700 km s -1 and a 
velocity resolution of 1 . 1 km s~ 1 . 3c279 was used as band- 
pass calibrator; MWC349, as flux calibrator (1.8 Jy); and 
0102+584 (1.1 Jy) and BLLAC (5Jy), as gain calibra- 
tors. Data calibration was performed within the MIR 
IDL package, while imaging and cleaning was done with 
MIRIAD. 

6 See http:/ /www. iram.es/IRAMES/mainWiki/Iram30mEfficiencies. 

7 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



Extremely broad RRL masers toward HW2 



3 



3. RESULTS 

In Figure 1 (upper panels), we present the full spectra 
measured toward Cepheus A HW2 at the wavelengths of 
the H40a, H34a and H31o: RRLs, and smoothed to a ve- 
locity resolution of ~6kms _1 . Since the HW2 radio jet 
is unresolved in the PdBI and SMA images (size of <2"; 
Rodriguez et al. 1994; Curicl ct al. 2006), the H40a and 
H31a spectra were averaged within the 5" -beam of the 
PdBI and SMA observations. The measured RRL spec- 
tra (Figure 1) show three different features: i) a strong 
slope due to the increase of the dust continuum emission 
with frequency (S v oc i^ Q ); ii) a forest of narrow molecular 
lines arising mainly from the HC source in HW2 (Martin- 
Pintado et al. 2005; Jimenez-Serra et al. 2009); and iii) 
three faint and extremely broad features extending in ve- 
locity from ^-500 km s -1 to ^600 km s -1 (zero- intensity 
lincwidths of ^1100 km s -1 ). 

From the slope in the RRL spectra, we derive spec- 
tral indeces for the continuum emission of ~2, consistent 
with that calculated by Comito et al. (2007). The PdBI 
and SMA indeces were obtained after subtracting the 
SEDs of the bandpass calibrators 3C454.3, 3C273 and 
3C279. The derived continuum (dust+free-free) level at 
the frequencies of the H40ai and H31a RRLs is ^0.12 
and 0.69 Jy, respectively. For the H34a line, this level 
was affected by an anomalous diffraction within the El 
receiver, although this did not affect the observed contin- 
uum slope and molecular line flux. The 2 mm continuum 
level (~0.37Jy) is estimated from the measured 3 mm 
and 1.4mm continuum fluxes (Column9 in Table 1). 

From Figure 1 , it is clear that molecular line confu- 
sion becomes an issue in this kind of studies, since mul- 
tiple line blending could generate broad features mimic- 
ing broad emission from RRLs. Molecular lines from 
species such as S0 2 , C 2 H 5 CN, HC 3 N or C 2 H 5 OH are in- 
dcntificd at frequencies close to the RRLs (lower panels 
in Figure 1). However, their linewidths are narrow (~4- 
8kms _1 ), and the resulting blending features do not ex- 
ceed ^lOOkrns -1 . Note that other blending features are 
also found across the RRLs spectral bands, with similar 
widths. It is then unlikely that line blending is resposi- 
ble for the extremely broad emission, with zero-intensity 
linewidths of ^1100 km s" 1 , detected toward Cepheus A 
HW2. The probability to detect three similar (and ex- 
tremely broad) features as a result of line confusion at 
three different frequency ranges, is also very small. We 
thus conclude that these features are likely associated 
with the emission of the H40a, H34a and H31a RRLs, 
formed in the high- velocity HW2 ionized jet (Section 5). 

Figure 1 (lower panels) shows a zoom-in of the 
continuum-subtracted RRL spectra, smoothed to a ve- 
locity resolution of 15-17 km s . The continuum emis- 
sion was subtracted by fitting a polynomial function of 
order 1. 

Although less clear for the H40a line, the extremely 
broad RRL profiles are asymmetric and have a similar 
kinematical structure with two broad (red- and blue- 
shifted) velocity components. The two-component Gaus- 
sian fits of this emission are shown in Figure 1 (see solid 
and dashed lines). The derived parameters for the red- 
(R) and blue-shifted (B) components are given in Table 1. 

Table 1 also reports the contribution from only the 
free-free continuum emission at 3, 1.9 and 1.4 mm 



(Column 10), derived by extrapolating the HW2 fluxes at 
centimeter wavelengths and by assuming an spectral in- 
dex ^0.6 (Rodriguez et al. 1994). These values, however, 
should be considered as upper limits since the actual free- 
free emission at millimeter wavelengths could be smaller 
than reported in Table 1 (Comito et al. 2007). If overesti- 
mated, these fluxes would imply an even stronger maser 
amplification effect for the observed RRLs (Section 4). 

The derived integrated intensities of the RRLs lie 
above the 9a Area level (Column 8), with a Area the in- 
tegrated intensity error in the spectra (see caption in 
Table 1). The derived peak and integrated intensities not 
only show a systematic trend to increase for increasing 
frequency (or decreasing quantum number n) for both ve- 
locity components, but to be brighter for the blue-shifted 
emission compared to the red-shifted component. This 
asymmetry in the RRL profiles can be explained as RRL 
maser emission in the HW2 jet (Section 5). 

Ignoring the results from the H40a blue-shifted com- 
ponent (its peak velocity and linewidth are factors of ~4 
and 2 larger than for H34a and H31a; Table 1), the av- 
eraged peak velocities and lincwidths are, respectively, 
— 50kms _1 and 360 km s -1 for the blue-shifted emission, 
and 360 km s -1 and 200 km s" 1 for the red-shifted gas. 

4. LINE-TO-CONTINUUM FLUX RATIOS: RADIO 
RECOMBINATION LINE MASERS IN HW2 

Martin-Pintado et al. (1989) showed that the genera- 
tion of RRL emission at millimeter wavelengths can suffer 
from non-LTE effects that lead to the formation of RRL 
masers in ionized winds. Since Cepheus A HW2 shows 
similar properties to the RRL maser object MWC349A, 
it is possible that the RRLs toward HW2 also form undcr 
non-LTE conditions. 

In Section 3, we have shown that there is a trend for 
the RRL peak and integrated intensities to increase with 
frequency. This behavior is expected for RRL emis- 
sion, since the integrated line-to-continuum flux ratios 
for RRLs (ILTRs, defined as AvT L /T c ; Martin-Pintado 
et al. 1989) increase with frequency as v 1A for opti- 
cally thin emission and LTE conditions. From Table 2 
(Columns 2 and 3), however, we find that the measured 
ILTRs for these lines clearly deviate from the LTE pre- 
dictions. The RRLs detected toward HW2 are therefore 
affected by maser effects. Like MWC349A (Column 4 in 
Table 2), the deviations from LTE for the ILTR of the 
H40a line are less prominent than those of H34a and 
H31a at higher frequencies. This is explained by the 
fact that the (3 n factors (the (3 n factor is proportional to 
the effective absorption coefficient of the RRL; Walmsley 
1990) are expected to be more negative for RRLs with 
n<35 than for those with n=40, at the electron densi- 
ties of ~l-5xl0 4 cm~ 3 derived toward HW2 (Curiel et 
al. 2006). 

5. MODELLING OF THE EXTREMELY BROAD RADIO 
RECOMBINATION LINE EMISSION 

We try to reproduce the observed RRL profiles to- 
ward HW2 by using the 3D radiative transfer model of 
Martin-Pintado et al. (1993, 2002, 2011), which includes 
the LTE departure coefficients, b n and /3 n , calculated by 
Walmsley (1990). The model considers an isothermal 
(~ 10 4 K) collimated, bi-conical ionized radio jet with a 
semi-opening angle of 18° and an inclination angle with 
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Table 1. Derived parameters of the H40o, H34a and H31a RRLs to- 
ward Cepheus A HW2. 



RRL 


v (MHz) 


A (mm) 




Tl (Jy) 


vlsr (kms 1 ) 


Ad (kms 1 ) 


T L Av (Jykms" 1 ) 


Sl ot (Jy) a 


S f J (Jy) b 


H40a 


99022.96 


3.0 


B 


0.0029 (0.0010) c 


~-210 


~620 


1.9 (0.1) d 


0.12 


0.059 








R 


0.0034 (0.0010) 


~368 


~190 


0.68 (0.05) 






H34a 


160211.52 


1.9 


B 


0.033 (0.009) 


60 


~410 


14.4 (0.7) 


0.37 c 


0.082 








R 


0.024 (0.009) 


~400 


~180 


4.4 (0.5) 






H31q 


210501.77 


1.4 


B 


0.055 (0.006) 


36 


~313 


18.5 (0.4) 


0.69 


0.099 








R 


0.035 (0.006) 


-316 


~240 


8.9 (0.4) 







a Total continuum (dust + free-free) emission flux measured toward HW2. 

b Free-free continuum flux obtained by assuming a spectral index of ~0.6, as derived from the VLA data of Rodriguez et al. (1994). 

c The error in the RRL peak intensity corresponds to the lcr noise level in the spectra (lmjy, 9mjy and 6mjy for the H40a, H34o and 
H31a RRLs respectively; see lower panels in Figure 1). 

d The error in the integrated intensity flux of the RRLs is calculated as <TA rfLa =laXy/&v X Av, with 8v the velocity resolution of ~15- 
17kms _1 in the RRL spectra (see Section3), and Av the linewidth derived in the Gaussian fit (Column 7 in this Table). 

c Derived from the continuum flux measured with the PdBI and SMA at 3mm and 1.4mm . 




500 
v lsr ( kms ) 

Fig. 2. — Upper panel: Morphology of the 3.6 cm radiocontinuum 
emission predicted by our model for the HW2 radio jet, smoothed 
to the angular resolution (0.25" xO. 18") of the VLA observations 
of Rodriguez et al. (1994). The X and Y axis show the offsets 
(in arcseconds) with respect to the major and minor axis of the 
jet. Contours are the same as in Rodriguez et al. (1994), i.e. - 
3, 3, 6, 9, 12, 15, 20, 30 and 50 times the rms (50/xJy) in the 
VLA 3.6 cm image. Beam is shown at the lower left corner. Lower 
panels: Spectra of the H40a, H34a and H31o lines predicted by 
our model assuming LTE (dashed lines) and non-LTE conditions 
(solid lines), overlaid on the RRL profiles measured toward HW2 
(black histograms). Numbers indicate the factors by which the 
intensities of the predicted spectra have been multiplied in order 
to qualitatively compare the observations with the modeling (see 
Section 5). 



Table 2. Integrated Line-to-Continuum flux ratios (ILTRs). 



RRL 


HW2 a (kms" 1 ) 


LTE b (kms" 1 ) 


MWC349 (kms" 1 ) 


H40a 


43 


~30 


39 c 


H34a 


229 


~43 


110 d 


H31a 


280 


-58 


215 c 



a Derived from the results shown in Columns 8 and 10 of Table 1. 
Calculated for optically thin continuum emission, T* = 10 4 K and 
iV(He+)/7V(H+)=0.08. 

c ILTRs of the H41a and H31a RRLs measured toward MWC349A 
by Martin-Pintado et al. (1989). 

d ILTRs of the H34a RRL obtained toward MWC349A by Thum, 
Martin-Pintado & Bachiller (1992). 



respect to the line-of-sight of 52°. This inclination an- 
gle is similar to that derived by Patel et al. (2005, of 
62°±10°) or Vlemmings et al. (2010, of 56°). However, 
the semi-opening angle (18°) is a factor of >2 larger than 
that derived by Rodriguez et al. (1994, of -7.5°). The 
determination of this angle is subject to large uncertain- 
ties since the radiocontinuum emission of the HW2 jet is 
variable, and its semi-minor axis is unresolved (Curiel et 
al. 2006). Despite the 18° semi-opening angle, our model 
reproduces well the radiocontinuum fluxes at centimeter 
wavelengths, and the extremely broad RRL profiles mea- 
sured toward HW2 (see below). 

The radiocontinuum spectrum of HW2 is fitted with a 
mass loss rate of 3.2xlO _6 M0 yr _1 . Although this rate 
is a factor of 4 larger than that estimated by Rodriguez et 
al. (1994, ^8xlO _7 M Q yr _1 ), the predicted radiocontin- 
uum fluxes at centimeter wavelengths agree well (within 
13%) with those reported by these authors (Table 3). In 
Figure 2 (upper panel), we show the morphology of the 
3.6 cm (8.44 GHz) radiocontinuum emission predicted by 
our model for the HW2 radio jet. This image is very 
similar to that observed with the VLA (Rodriguez et al. 
1994). The derived free-free continuum fluxes at 3, 1.9 
and 1.4mm (0.058, 0.079 and 0.094 Jy, respectively) are 
also consistent with those shown in Table 1. 

The assumed electron density distribution in the HW2 
jet varies as r -211 , with r the distance to the central 
protostar extending out to 0.41" (285 AU at a distance 
of 700 pc; Reid et al. 2009). The electron density at the 
inner radius of 6.7 AU is 2.7xl0 8 cm -3 , sufficient to ob- 
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Table 3. Comparison between the observed and predicted radiocon 
tinuum fluxes toward cepheus a hw2. 



Frequency 
(GHz) 


Observed S„ a 
(mjy) 


Predicted S„ h 
(mjy) 


Relative error 

(%) 


1.49 


3.4 


3.3 


1.6% 


4.86 


7.5 


7.6 


1.1% 


8.44 


9.8 


11.1 


12.9% 


14.9 


15.8 


16.3 


3.2% 


43.3 


35 


33.4 


4.5% 



a From Rodriguez ct al. (1994). 
b Predicted by our model 



for 



rate of 



3.2x10 Mq yr , and a semi-opening angle of 18°. 



tain the maser effect in the RRLs. 

To explain the extremely broad linewidths and the lack 
of emission at systemic velocities, we have considered 
that the ionized gas in the jet is accelerated constantly 
to reach a terminal velocity of 500 km s -1 (Curiel et al. 
2006) at 35 AU from the protostar. The model also in- 
cludes electron impact (pressure) broadening and a tur- 
bulent velocity of 60kms~ 1 . However, pressure broad- 
ening is expected to be negligible for the observed RRLs 
(Keto et al. 2008). 

Figure 2 (lower panels) shows the RRL profiles pre- 
dicted under LTE (dashed lines) and non-LTE conditions 
(solid lines), overlaid on the observed RRLs. To match 
the observations, we needed to (red-)shift the predicted 
RRLs by 160 km s -1 with respect to the systemic velocity 
of the cloud (-lOkms" 1 ; Martin-Pintado et al. 2005). 
Water masers toward HW2 also show significant shifts 
in their peak velocities, associated with the dynamics 
of the HW2 rotating disk and wide-angle outflow (Tor- 
refies et al. 1996, 2011). However, their velocity spans 
(~40kms _1 ) are too small to explain the 160 km s -1 
shift. The morphology and kinematics of the inner re- 
gions in the jet are unknown. Asymmetries in the den- 
sity, temperature and kinematics between the red and 
blue lobes of the jet could give rise, due to the maser 
phenomenon, to large asymmetries in the RRL profiles, 
as expected from the larger spatial scales needed for 
RRL maser amplification compared to molecular masers 
(sec below) . Future modelling will explore whether these 
asymmetries are responsible for the 160 km s -1 shift ob- 
served in the RRLs. 

From Figure 2, we find that the model qualitatively 
reproduces the two-component line profiles, with zero- 
intensity linewidths of ^1000 km s -1 , of the RRLs to- 
ward HW2. It can be seen that maser amplification 
(solid lines) occurs at radial velocities smaller than 
±300 km s" 1 with respect to the systemic velocity of 
the jet. This is due to i) high-enough electron densi- 
ties at those velocities (3xl0 5 -4xl0 7 cm~ 3 ; Strelnitski, 
Ponomarev & Smith 1996); and ii) large enough coher- 
ent lengths, which lead to substantial amplification of 



the radiation along a velocity-coherent path (Ponomarev 
1994). For larger velocities (>±300kms _1 ), the electron 
densities are also high (3x 10 5 -8x 10 6 cm -3 ). However, 
the coherence length is very small making the maser ef- 
fect negligible (Ponomarev 1994). Maser amplification 
is maximum at radial velocities of ±140-200 km s _1 with 
respect to the systemic velocity of the jet, because the 
electron densities at these velocities correspond to the op- 
timum values for maser amplification (i.e. 6.8xl0 6 cm~ 3 
for H40a, 1.7xl0 7 cnr 3 for H34a, and 3.0xl0 7 cnT 3 
for H31a; Figure 8 in Strelnitski et al. 1996). Toward 
MWC349A, Martin-Pintado et al. (1994) also reported 
the detection of RRL maser emission at velocities very 
different from the ambient cloud velocity (±60kms~ 1 ), 
and predicted that they arise from ionized outflowing gas 
(Martin-Pintado et al. 2011). 

Although the model reproduces well the intensity 
of the H40a line, it fails to predict the intensities of 
the H34a and H31a lines (note that these intensities 
have been multiplied in Figure 2 by factors of 5 and 3, 
respectively). This problem of unmatched intensities for 
RRLs with n<40, was already noted by Martin-Pintado 
et al. (1989) for the MWC349A RRL maser object, and 
could be due either to small density and temperature 
inhomogeneities in the stellar wind, or to uncertainties in 
the (3 n coefficients (Mart in- Pint ado et al. 1993). Despite 
the unmatched intensities, our model reproduces the 
asymmetry observed in the RRLs (Section 3), which can 
only be explained by RRL maser effects. All this sug- 
gests that the extremely broad features detected in the 
millimeter spectra toward Ccpheus A HW2, are associ- 
ated with RRL masers generated in the HW2 ionized jet. 

In summary, we report the detection of the H40a, 
H34a and H31a RRLs toward the high-velocity ion- 
ized jet in Cepheus A HW2. This emission shows ex- 
tremely broad line profiles with zero-intensity linewidths 
of ^llOOkms -1 . The derived ILTRs significantly de- 
viate from those in LTE, suggesting that these lines 
are RRL masers formed in the HW2 ionized jet with 
a expanding velocity of 500 km s -1 . Together with 
MWC349A, the discovery of RRL masers in HW2 sug- 
gests that this mechanism could be a common feature in 
dense UC HII regions (Martin-Pintado 2002). 
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